Vortex dynamics and second magnetization peak in PrFeAsO0.60F0.12 superconductor 
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We have studied the vortex dynamics in the PrFeAsO0.60F0.12 superconducting sample by dc 
magnetization and dynamic magnetization-relaxation rate (Q) measurements. The field dependence 

f/-\ of the superconducting irreversible magnetization M s reveals a second magnetization peak or fishtail 

,_^ effect. The large value of Q is an indication of moderate vortex motion and relatively weak pinning 

energy. Data analysis based on the generalized inversion scheme suggests that the vortex dynamics 

£SJ can be described by the collective pinning model. The temperature dependence of the critical 

current is consistent with the pinning due to the spatial variation in the mean free path near a 
lattice defect (81 pinning). The temperature and field dependence of Q indicates a crossover from 

^H elastic to plastic vortex creep with increasing temperature and magnetic field. Finally, we have 

constructed the vortex phase diagram based on the present data. 
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I. INTRODUCTION 



The relatively high superconducting transition temperature T c up to 55 K and phenomenally high upper critical 
O . field H C 2 along with lower anisotropy and larger coherence length compared to the high-T c cuprates have made 
the pnictides interesting candidates for possible technological applications [H [3J. Another important feature of 
^ these compounds is the presence of a peak in the field dependence of critical current density, which is associated 
with a fishtail or second magnetization peak (SMP) effect in the M(H) curve, as observed in several cuprates and 
conventional superconductors [3TU2J. The appearance of SMP in the M(H) curve is not a universal phenomenon in 
pnictides. Whether in a particular material the SMP would appear is very much sample specific [T3HT6] . Most of the 
*^ theoretical approaches agree that the temperature dependence of the field H p , at which the SMP has its maximum, 
is related to the nature of vortex pinning and corresponds to a crossover between two different regimes of vortex 
lattice. Even though there are a lot of studies dealing with this phenomenon, the mechanism and origin of this effect 
are still much debated partly because the proposed models are system specific. 



7— \ So far, the studies of vortex dynamics have been mainly focused onto the doped 122 family AFe2As2 (A = Ca, 

Ba, Eu, etc.) as sizeable single crystals are available in this group of pnictides. In the 122 compounds such as 
BaFe 2 _ 3; Co a ;As2 [Ml [14], Cai_ a; Na a ;Fe2As2 [15] . and Bai_ a; K ;r Fe2As2 [IS] the vortex dynamics has been described 
by the collective pinning model. In these compounds, the SMP in the M(H) curve signifies a crossover from elastic 
to plastic vortex creep regime. In contrast, studies of Bitter decoration, small-angle neutron scattering and magnetic 
force microscopy in 122 compounds have revealed highly disordered vortex-glass phase with a short-range hexagonal 
order [T7H19] . Furthermore, these studies suggest that the vortices remain in the single- vortex pinning limit even 
at high-magnetic fields up to 9 T. Kopeliansky et al. 20 proposed that the SMP in Ba(Feo. 925000.075)2^82 is 
associated with a vortex structural phase transition from rhombic to square lattice occuring at field and temperature 
• • corresponding to the minimum point of the magnetic-relaxation rate. In LiFeAs too, the SMP has been attributed to 
a vortex structural phase transition [21]. In BaFei.82Nio.isAs2, though a fishtail effect appears in the M(H) curve, 
it is not associated with the crossover in a vortex pinning regime within the collective pinning scenario |22j . On the 
other hand, in Ca(Fei_ 2 ;Co x )2As2 the SMP is absent and the analysis of temperature- and field-dependent magnetic 
relaxation data suggests that the vortex dynamics is consistent with the plastic creeping model [23 . In contrast to 
the above, there are only few studies on vortex dynamics in the 1111 family of i?FeAs(0,F) (R = La, Ce, Pr, Nd, Sm 
etc.) compounds. Investigation of vortex dynamics by magnetization-relaxation measurements in SmFeAsOo.gF .i 
[2~4] and NdFeAsOo.gFo.i [55] [55] suggests that the vortex behavior is consistent with weak collective pinning and 
creep. Local measurements by magneto-optical imaging and microscopic Hall sensors in underdoped PrFeAsOo.g 
and NdFeAsOo.gFo.i indicate that the pinning is due to the collective creep and the origin of the SMP is due to a 
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order-disorder transition of the vortex lattice [27J . 



As the vortex dynamics of the pnictide superconductors is sensitive to doping, anisotropy, etc. one has to study 
each system separately as a function of these parameters in order to understand the origin of the SMP. In this report, 
we have investigated the vortex dynamics and origin of SMP in the optimally-doped PrFeAsO0.60F0.12 sample by 
dc magnetization measurements, and temperature and field dependence of dynamic magnetization-relaxation rate 
Q. Data have been analyzed by the vortex collective pinning model and the method of generalized inversion scheme 
(GIS) [251 US]- Analysis of the temperature and field dependence of magnetic relaxation suggests that the vortex 
dynamics can be described by the collective pinning and creep, and a crossover from the elastic to plastic creep 
regime with increasing temperature and magnetic field. 



II. EXPERIMENTAL TECHNIQUE 



The polycrystalline sample used in this study was prepared by solid-state reaction method. The details of the sample 
preparation, powder x-ray diffraction analysis and magnetic characterization have been discussed in our earlier reports 
[301 131) . The phase purity of the sample was examined from x-ray and energy dispersive x-ray analysis [301 131] . Both 
the zero field cooled (ZFC) and field cooled susceptibilities start to deviate from the normal behavior below 48 K due 
to the appearance of a diamagnetic signal, which is close to the zero-resistance temperature (30]. At 4 K, the shielding 
and Meissner fractions are calculated to be 88% and 35%, respectively, from the dc susceptibility data at H — 10 Oe. 
The shielding fraction of the sample is ~ 100% at 4 K determined from the real part of the ac susceptibility curve 
at H ac — 3 Oe. The dynamic magnetization-relaxation measurements were carried out in a VSM-PPMS (Quantum 
Design) using the following protocol. The sample was cooled down to the desired temperature in the ZFC mode and 
then the magnetic field is swept and we measured the magnetic moment following the routes: — > H max — > — 1 T — > 
0, with different field sweeping rates dH/dt = 40, 80 and 160 Oe/s. The magnetization-relaxation rate Q is defined as 



Q 



d\nj s 



dln(AM) 



dln(dH/dt) d\n(dH/dt)'' 



(1) 



where j s is the transient superconducting current density and AM 
magnetization for dH/dt <0 (dH/dt > 0). 



(M+ - M_), M+(M-) is the branch of the 




FIG. 1: (color online) (a) Magnetization hysteresis M(H) loop of the PrFeAsO0.60F0.12 sample at different temperatures. Inset: 
the field-decreasing branch of hysteresis loops at 2, 3 and 5 K showing the effect of flux jump in the form of small kinks due to the 
thermomagnetic instability, (b) M(H) loop at 53 K, just above the superconducting transition temperature. Arrows indicate 
the direction of increasing and decreasing magnetic field, (c) The superconducting irreversible signal M S (H) at 22 K exhibiting 
the second magnetization peak. M S (H) curve is obtained by subtracting the paramagnetic component M p = (M+ + M_)/2, 
where M+(M-) is the branch of the magnetization for dH/dt < 0(dH/dt > 0), from the measured M(H) loop. 



III. RESULTS AND DISCUSSION 



A. Magnetization-hysteresis loop and dynamic magnetic relaxation 
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FIG. 2: (color online) Field dependence of AM at different temperatures with field sweeping rate 160 Oe/s (solid line) and 40 
Oe/s (dotted line) and the dynamic magnetization-relaxation rate Q (red line). Left axis shows the magnetization and right 
axis is for the dynamic magnetization- relaxation rate. Arrows indicate the onset field, H on , and the peak field, H p , of the 
second magnetization peak. 

Figure fT] shows the magnetization hysteresis loops M(H) of the PrFeAsO0.60F0.12 sample in the superconducting 
and normal state (just above T c ). One can clearly see that the superconducting hysteresis loop arises from the flux 
gradient produced by the pinning of flux lines. A small section of the field-decreasing branch of hysteresis loops at 
2, 3 and 5 K is shown in the inset of Fig. flVa), wherein the signature of flux jumps is visible in the form of small 
kinks. Usually, flux jumps in superconductors occur in the low-field and low-temperature regime due to the effect 
of thermomagnetic instability (TMI) [55] • Due to the thicker sample used in our study, the effect of flux jump is 
not strong enough though it is detectable as shown in the inset of Fig. Ilia). The observed M(H) curve in the 
superconducting state is the sum of a superconducting irreversible signal (M s ) and a paramagnetic component (M p ). 
Analysis of dc susceptibility, x(T), in the normal state shows that the paramagnetic component of magnetization 
comes from the Pr 3+ magnetic moments 13 f j . As the sample has a large magnetic background in the normal state, 
we have measured magnetic hysteresis loop at 53 K, just above T c as shown in Fig. [lib). Figure njb) shows no 
observable hysteresis, thereby excluding the presence of irreversible magnetization in the sample due to the free 
iron nanoparticle or magnetite traces. Figure ITVc) shows that the superconducting irreversible signal M S (H) at 22 
K exhibits a second magnetization peak. AI^yH) curve is obtained by subtracting the paramagnetic component 
M p = (M + + M_)/2 from the measured M + and M_ branches of M(H) loop. In a polycrystalline sample, the gap 
AM in the magnetization loop can be split into intergranular (global) and intragranular (local) parts |33j . In the 
low-magnetic field region, AM is predominantly caused by the intergranular current, but in the high-field region AM 
results largely due to the intragranular current. This has been confirmed from the magnetization loop measurements 



on the bulk and powder sample in cuprates and pnictides [33J [35] . From the study of remnant magnetization as a 
function of maximum applied field in Sm and Nd iron oxypnictides, Yamamoto et al. [36] showed that the intergrain 
current persists only up to a few hundred Oe. Since our study goes beyond this low-field region, we may assume that 
AM originates mainly from the intragrain current density. 

In order to investigate the mechanism of vortex pinning and origin of the SMP in PrFeAsO0.6F0.12 sample, we 
have followed the dynamic magnetic-relaxation technique. Figure [2] shows the field dependence of AM measured 
with the field sweeping rates of 160 and 40 Oe/s, and the value of Q at different temperatures was calculated using 
equation (HI. The arrows in the figure point to the characteristic fields: the onset field, H on , and the peak field, H p , 
of SMP. The difference between AM measured at 160 and 40 Oe/s can easily be distinguished, which indicates a 
giant vortex creep as observed in the cuprate superconductors. It is also clear from the figure that both H on and H p 
are dependent on the magnetic field sweeping rate and their values decrease with the decrease in field sweeping rate. 
Similar vortex creep has also been observed in single crystals of Ba(Fei_ K Co a )2As2 d3J[H] and Bai_ :E K 3 ;Fe2As2 [IS], 
and polycrystalline SmFeAsOo.9F .i sample [24]. The value of Q at 1 T and 7 K is as large as 6% which is comparable 
to that observed in cuprate superconductors (e.g., 4% in YP^CusOy [37]) but one order of magnitude larger than 
that for MgB2 [35]. There are several important features that can be identified from the field dependence of Q: 
(1) H on is weakly temperature dependent in the high-T region but increases rapidly with decreasing temperature 
exceeding the limit of our measuring range (8 T) below 5 K. (2) In the low-temperature region (left panel of Fig. 
[2]), a peak in Q{H) is observed where the field dependence of AM changes its slope rapidly and the peak position 
shifts to lower field with increasing temperature. The low-field peak in the Q(H) curve has also been observed 
in Ba(Fei_a;Coj;)2As2 single crystal and has been explained as the signature of a crossover between two different 
regimes of vortex dynamics 13j. (3) In the high- temperature region (right panel of Fig. pi), there is a minimum in 
the Q(H) curve which falls in between the fields H on and H p . However, in case of Ba(Fei_ :E Co :r )2As2 single crystal 
this minimum roughly matches with the second peak position, H p , in the J{H) curve [13j . 

Figure [3] shows the temperature dependence of Q for different magnetic fields. Q shows a bell-like or bump-like 
shape in the intermediate-temperature region for low field. With increasing magnetic field, the bump-like feature 
shifts toward lower temperature and disappears above 4 T. For all magnetic fields, Q(T) passes through a minimum 
at a temperature T m as indicated by arrows in Fig. [3] At 0.5 T, T m « 30 K (~ 0.63 T c ) and with increase of magnetic 
field T m shifts toward lower temperature. Similar bump-like shape in low field has also been observed in the Q{T) data 
of Ba(Feo.92Coo.o8)2As2 single crystal [T3J and in the remanent magnetization relaxation data of Ba(Fe .g3Coo.o7)2As2 
[H] and Cao.25Nao.75Fe2As2 [IS] single crystals. The bump-like shape disappears and a plateau appears in Q(T) or 
S(T) for H above 1 T and 3 T for Cao.25Nao.75Fe2As2 [H] and Ba(Feo.92Coo.os)2As2 [13], respectively. In all these 
cases, this behavior has been attributed to the effect of collective pinning of vortices. 




FIG. 3: (color online) Temperature dependence of the dynamic magnetization-relaxation rate Q in different magnetic fields. 
Arrows indicate the temperature T m at which Q(T) shows a minimum. 



B. Analysis based on the vortex collective pinning model 
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FIG. 4: (color online) (a) Temperature dependence of the T/Q in different magnetic fields. The elastic to plastic crossover 
temperature, T ep , was determined from the peak position of the T/Q vs T curve for different magnetic fields. Arrows indicate 
(b) -dlnj B /dT vs Q/T at 7 K in different magnetic fields. The best fit line (red) has slope C — 21.62±1.2. 



T ep for 2 and 4 T 



In a type-II superconductor, the vortices normally move over the effective pinning barrier U(j s ,T,H e ) by thermal 

activation with an average velocity v = Vq exp I ^' ^ "' J , where Vq is the attempt hopping velocity and H e is the 

actual local magnetic induction. The electric field, E, induced by this vortex motion is given by [39] . 



E = vqH exp — 



U( Js,T,H e ) 



(2) 



In a system with randomly distributed weak pinning centers, the current dependent effective pinning barrier can be 
written in a general form [ID] 



U(j s ,T,H e ) 



U c (T,H e ) 



Jc (T,H e ) Y (TH 
3s(T,H e ) 



(3) 



where /x, U c , and j c are the glassy exponent, intrinsic characteristic pinning energy, and the unrelaxed critical current 
density, respectively. According to the collective pinning theory, for a three-dimensional lattice the values of exponent 
\i are 1/7, 3/2 and 7/9 for the single vortex, small bundles and large bundles of vortex motion, respectively [41] . 
For fi = — 1, equation pi describes the Kim-Anderson model [22] and in the limit /i — >0, it reduces to the Zeldov 
logarithmic model [42 . Any value of /i, positive, negative or zero is physically meaningful. A negative /i corresponds 
to a finite dissipation in the small current limit and plastic vortex motion, while a positive \i corresponds to a vanishing 
dissipation in the small current limit and elastic vortex motion |43) . From the general equations (l)-(p|, Wen et al 
[2"9"] derived the following equation: 



T 



Q(T,H e ) 



U c {T,H e 

kn 



li(T, H e )CT 



(4) 



where C = In [2voH e /l(dH e /dt)] is a parameter that is weakly temperature dependent and I is the lateral dimension 
of the sample. Figure [i|a) shows the T/Q vs T curves at different magnetic fields. The T/Q curves go through a 
peak, the position of which varies with the applied field. The left side of the peak denotes a positive slope of the 
curve which corresponds to a positive /i [equation H — provided U C (T) has no strong temperature dependence — and 
thereby elastic vortex motion while the right side of the peak has a negative slope which suggests a plastic vortex 
motion. Thus, as the strength of the applied field is increased, the low-temperature region over which the elastic 
vortex motion prevails is shortened. By extrapolating the curve T/Q down to zero temperature, we have deduced 
the value of £/ c (0) at 0.5 T from Fig. |4|a). The obtained value of U c (0)/kB is about 36 K, which is comparable to 



those for SmFeAsOo. 9 F .i (C7 C (0) ~ 40 K) [23] and Ba(Fei_ a; Co a: )2As2 (t/ c (0) - 98 K) [13]. However, these values 
are order of magnitude smaller compared to U c (0) ~ 300 K for YBCO thin films [29] and beyond 3000 K for MgB2 
[33], implying a quite small characteristic pinning energy for the pnictide superconductors. Assuming U C (T) is not a 
strong temperature-dependent function, the value of \iC can be obtained from the slope of the T/Q vs T curve. At 
H = 0.5 T, fiC is evaluated to be 11.57±0.64. To estimate C, one can make use of the equation [25] 



dT 



dAnjc 
dT 



C 



(5) 



which is valid in the temperature region where j s is not affected by the quantum creep at low temperatures and lnj s 
follows an almost linear decrease with T. —d\nj s /dT vs. Q(T)/T has been plotted at 7 K for five different magnetic 
fields as shown in Fig. pjb). The slope of the best-fit straight line provides the value of C = 21.62 ± 1.2. The 
parameter C can also be evaluated from the equation [2SJ 



C = lim 



• 1 d In j s 



t->o Q d\nT 



(6) 



The extrapolation of (— l/Q)(dlnj s /dhiT) vs. T curve to K gives C = 21.8. The values of C obtained 
from these two procedures are close to each other and slightly smaller than that reported for Ba(Fei_ a; Co :r )2As2 
single crystals (C= 28.56) [T3J but close to that for YBa2Cu40g (C= 19) [25]. From the values of \iC and C, 
we have estimated the value of /x to be 0.53. The positive value of (i indicates an elastic vortex motion. In the 
later section, we will show that this value of fi is consistent with that estimated from the generalized inversion scheme. 



C. Analysis based on generalized inversion scheme 
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FIG. 5: (Color online)(a) Temperature and field dependence of the unrelaxed AM c (T,H e ) evaluated by following GIS. (b) 
Temperature dependence of jc(T)/j c (0) in different magnetic fields along with the theoretical prediction for 51 and 8T C pinning. 



In order to evaluate the temperature and field dependence of the unrelaxed critical current density j c (T,H e ) and 
the corresponding characteristic pinning energy U c (j s ,T, H e ) directly from relaxation data, we have followed the 
generalized inversion scheme proposed by Schnack et al [28 and Wen et al [29 . Though the present sample is 
polycrystalline in nature, scanning electron microscope images reveal large, oriented, and well-connected plate-like 
crystallites. The average grain size in this sample is larger than those normally reported in oxypnictide polycrystalline 
samples and the plate-like grains have a tendency to orient their a6-plane along the broad face of the sample [3T] . All 
these led us to venture for the GIS analysis on the relaxation data of our polycrystalline sample. The GIS scheme is 
more general as it does not require a priori assumption about the explicit temperature or field dependence of U. The 
basic postulates of GIS are: (1) thermally activated flux motion can be described by equation p); (2) the activation 
energy can be expressed as product of two general functions g(t = T/T c ,H e ) and f\j s (T,He)/j e (T,H e ),H e ], i.e., 
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FIG. 6: (color online) Current dependence of the effective pinning energy U(T,H e ) in different magnetic fields evaluated by 
following the method of GIS. The solid line is a theoretical fit of equation (13j) to the U(j B ,T) data for H = 0.5 T in the 
temperature region 2-20 K. Using the obtained best fit parameters U c (0) = 48 K and [i = 0.59 this curve has been extended 
for T > 20 K. 



U(j s , T, H e ) = U c {0, H e ) x g(t, H e ) x / [ Js (T, H e )/ 3c (T, H e ), H e ] with /(l, H e ) = 0, since U( 3c , T, H e ) = by definition 
of the critical current density j c ; and (3) g(t, H e ) oc \j c (T, H e )/j c (0, H e )] p G(t), here G(t) and p depend on the specific 
pinning models [29] . According to GIS, the temperature dependence of the true critical current density j c at a 
particular field can be evaluated from the following integral [251 [5S] 



jc{T) = Jc(0) x exp 



CQ(T')[1 



d In G(T') i 
dlnT' 1 



dlnj B (T') 
dlnT' 



l+pQ{T')C 



dT 



(7) 



where j c (0) is the true critical current density at K. In the above equation, Q(T) and j s (T) are the measured values. 
In order to apply the above procedure, one must know the value of p as well as the function G(t). For a two-dimensional 
pancake system, p = 1 and G(t) = y/{\ + t 2 )/(l — t 2 ) with t = T/T c . Similarly, for a three-dimensional single vortex, 
p = 0.5 and G{t) = (1 + i 2 ) 5 / 4 /(l - t 2 ) 1 / 4 (Ref.[25]). In the underdoped PrFeAsO . 9 , it has been shown that in the 
low-field region the vortex pinning is in the three-dimensional single- vortex limit [27] . Therefore, we have taken p = 
0.5 and C ~ 22 to calculate the values of j c (T) following the method of GIS. As the supercurrent density j s ex AM, 
therefore, we have plotted AM C (T) in Fig. ISFa) for different magnetic fields. In type II superconductors, the pinning 
originates from two basic mechanisms: due to the spatial fluctuation of superconducting transition temperature T c 
(<5T c -pinning) and due to the spatial variation in the mean free path I near a lattice defect (iM-pinning). In the 
single vortex regime, the theoretical curve for Spinning is given by j c (T)/j c (0) = (l + i 2 ) _1 / 2 (l — £ 2 ) 5 / 2 and that 
for the £T c -pinning is given by j c (T)/j c (0)=(l - t 2 ) 7 / 6 (l + t 2 ) 5 / 6 (Ref. [1511351). In Fig. [HJb), we have shown the 
temperature dependence of experimentally derived j c (T)/j c (0) values together with the theoretical predictions for 
two basic pinning mechanisms. It is clear from the figure that the experimentally determined values are close to the 
^/-pinning mechanism. This is also consistent with the prediction that the flux pinning in the charge-doped pnictides 
can be described by the mean- free path fluctuations introduced by the dopant atoms [37] • Figure [6] shows U(j s , T, H e ) 
obtained by GIS at three different fields for the PrFeAsO0.60F0.12 sample. From the fitting of equation (pi) to the 
data for 0.5 T, the parameter (7 C (0) is evaluated to be 48 K which is close to the value (~36 K) determined from the 
analysis using collective pinning model. From the fit, we also get fi = 0.59, which is close to that obtained (fi — 0.53) 
from the T/Q vs. T curve. It should be noted that the values of /i determined here reflect just an averaged one, 
which, in principle, is also current dependent. 



IV. PHASE DIAGRAM AND SECOND MAGNETIZATION PEAK 

In Fig. [7| we have plotted the temperature dependence of the characteristic fields, viz., H on , H p , H ep (the elastic 
to plastic crossover field), H rr and H&. The H ep (T) values are obtained from the T/Q vs. T curves for different H 
(FigEJa)). The H ep (T) curve is located in between the onset, H on and the peak, H p of the SMP. The characteristic 
fields H on and H p follow a concave-shaped decrease with the increase in temperature. The SMP has been observed in 
the conventional as well as high-temperature superconductors and various mechanisms have been proposed to explain 
the origin of the same. Foremost among them are (i) a change in the dynamics of the vortex lattice — a crossover 
from the single- vortex pinning regime to the bundle pinning regime [1], (ii) a change in the vortex creep mechanism 
— a crossover from the elastic to plastic creep [H], (iii) a first-order phase transition from an ordered "elastically 
pinned" low-field vortex phase, the so-called Bragg-glass [49] to a high-field disordered phase characterized by the 
presence of topological defects [50j H2] and (iv) a structural vortex lattice transition (from rhombohedral to square 
[52"] or from triangular to square lattice [53j) well within the Bragg glass regime. 

For the present case, the absence of a mirror image correlation between AM and Q rules out the possibility of 




FIG. 7: (color online) Vortex phase diagram of PrFeAsO0.60F0.12 sample describing the temperature dependence of the onset 
field, H on , the elastic to plastic creep crossover field, H ep , the peak position H p , the irreversibility field H irr and the upper 
critical field H C 2- Data for Hi rr and H C 2 are taken from Ref. |31| and |30| . respectively. The dashed line and dotted line 
describe the power law fit (1 — ^-)" to H p and Hi rr data with n = 1.6 for H p and 1.7 for Hi rr . 



a crossover from the single- vortex pinning regime to the bundle pinning regime as the origin of SMP [3]. The T 
dependence of T/Q curve for different H (Fig. El) has unambiguously established the occurrence of a crossover in 
the flux dynamics, from elastic to plastic creep with increasing temperature. However, this crossover line, H ep , does 

on or H p of the SMP; rather, it lies in between H on and H p . 



It may be mentioned that in YaBaCu07_a, H, 



<:p 



not coincide with any of the characteristic fields, viz.,_ff, 

Also, there is no anomalous feature in the M(H) curve at H ep . 

coincides with H p [HE], whereas in Ndi.ssCeo.isCuC^-a the creep crossover occurs well above H p [5J. In case of 

the crossover in the flux creep behavior occurring below H p (T), both H p (T) and Hi rr (T) are controlled by plastic 

pinning and their temperature dependence are expected to be of the form (1 — ^-) n [13j[l4]. For Ba(Fei_ :r Co :r )2As2 

(x = 0.08) single crystal, H p (T) and Hi rr (T) lines follow such kind of temperature dependence with n = 1.7 and 

1.52, respectively [13]; for x = 0.07, H p (T) goes with n = 1.5 [14 . In the present case, we observed n = 1.6 and 1.7 

for the H p (T) and Hi rr (T) lines respectively, in agreement with the data for Ba-122 crystal. 



It is well accepted that a crossover from elastic to plastic vortex creep may accompany a quasi-order-disorder (OD) 
transition of the vortex lattice. The vortex phase of the superconductors is determined by the competition between 
the energy of thermal fluctuations, the pinning energy generated by the quenched disorder, E p , and the elastic energy 
of the vortex system, E e i [5, 49, 54, 55 . E p and E e i are directly related to the superconducting parameter such 
as penetration depth A, coherence length £, pinning parameter and the anisotropy factor. At low T, where thermal 
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FIG. 8: (color online) (a) The low temperature variation of H on (T) for T/T c < 0.5 in double logarithmic scale. In the low-T 
range, H on (T) oc 1/T 1 ' 1 . The solid lines represent a linear fit in two temperature intervals to emphasize the presence of a 
inflection-like behavior located between f2 and 15 K. (b) The variation of C/(~ T/Q) with H at 2 K on a log-log scale. Solid 
line describes the power law fit U(H) oc H 0A for the H domain corresponding to H on at T/T c < 0.5. 



energy is small compared to E p and E e i, the OD transition is roughly described by the equality [5, 6 

E p (T,H)=E el (T,H). 



(8) 



For static condition where no current flows in the sample, OD transition field is practically independent of T in the 
low-T region, far below T c where the superconducting parameters vary weakly with temperature [SHE]. In global 
magnetic measurements, the transition is relatively wide, and the transition field lies somewhere in between H on 
and H p , though H on is usually considered as the transition field. In cuprates, the occurrence of SMP is commonly 
associated with this OD transition [SHE]- Unlike cuprates, in the present case H on exhibits a concave-shaped increase 
with the decrease of T, as observed for other characteristic fields such as, H ep and H p . I^-^Sr^CuO^ with x = 
0.126 and similar doping, exhibits a broad SMP with characteristics that are strongly temperature dependent down 
to low temperature [52] [5j5] . This behavior was explained by considering the softening of the vortex lattice associated 
with the square to rhombic vortex lattice transition as the source for the SMP [30] [52] . Alternatively, the upward 
curvature in the T dependence of H on and H p in the low-T region may be explained as a dynamic effect caused by the 
finite current j s induced in the specimen during standard dc magnetization experiments, which reduces the effective 
pinning energy U(j s ,T,H) [57]. In a crude approximation, in the dynamic condition, the appropriate energy balance 
relation should be 

U(j a ,T,H)<xE el (T,H). (9) 

With U(j s ,T,H) oc T and E e i(T,H) ex A^iT -1 / 2 (independent of j s ) one obtains a T~ 2 dependence for H on , 

H on (T) ex A- 4 T- 2 (10) 

at low temperature. In the high temperature range (0.5 < T/T c < 1), the dependence of superconducting parameters 
on T is important to determine the origin of pinning. In this region, for SI pinning, H on (T) is expected to show a 

1/9 

peak and follow the expression H on (0) [l — (T/T c ) 4 ~\ below the peak [5]. Indeed, we have observed a weak peak 

at T ~ 0.67T C which may be related to the <W-pinning of the vortices. However, due to limited data points in a 
narrow temperature range, fitting could not be done unambiguously. 

For the low-T range (T/T c < 0.5), H on (T) plotted in a double logarithmic scale is shown in Fig. [8la). In this 
region, H on varies as T which is smaller than the exponent predicted in equation (10). In FeSei-^Te^, such a 

low exponent was attributed to the increase of U with increasing H, specific for the collective creep regime [58) . From 
a plot of [7(~ T/Q) vs H at 2 K (Fig. ^[b)) we observed a H 0A dependence of U(H) in the H domain corresponding 



to H on and at T/T c < 0.5, indicating U(T,H) varies as TH 0A . Combining this with equation QQJ) , and neglecting the 
T dependence of A in the low-T limit, one obtains H on {T) oc T~ ln , as observed in Fig. pMa). The observed upturn 
in U(H) below 2 T (Fig. Isfb)) may be attributed to the the effect of TMI appearing in the low-T domain [55]. There 
is another striking similarity of our result with those for FeSei_^Te K [581 and La2- :E Sr :E Cu04 crystals |57j ; we find an 
inflection-like point in the H on vs T behavior between 12 and 15 K (Fig. pja)). The location of this inflection point 
is found to be in good agreement with the results obtained for X(T) in underdoped PrFeAsOi_j, single crystal [59] 
and following Miu el al. [57] [55] we associate this feature with the T dependence of the superfluid density (n s ex A -2 ) 
in the case of two-band superconductivity. 

V. CONCLUSION 

In summary, we have investigated the vortex dynamics of the optimally-doped PrFeAsO0.60F0.12 sample by dc 
magnetization and dynamic magnetization-relaxation measurements. The field dependence of the superconducting 
irreversible magnetization reveals a second magnetization peak. In low-field and in the intermediate-temperature 
region, Q(T) exhibits a bell-like or bump-like shape. In this region, data analysis based on the vortex collective 
pinning model and GIS suggests that the vortex dynamics can be better described by elastic vortex motion. Analysis 
of the temperature and field dependence of Q suggests a crossover of the vortex dynamics from the elastic to plastic 
creep regime with increasing temperature and magnetic field. The temperature dependence of the critical current 
density is consistent with the pinning due to the spatial variation in the mean free path near a lattice defect (51- 
pinning) . Analysis of present data persuades for an order-disorder like transition of the vortex phase in our sample 
followed by a crossover from elastic creep to plastic creep in the vortex motion. 

VI. ACKNOWLEDGEMENT 

The authors would like to thank Prof. S. Das, N. Khan, A. Midya and A. Paul for technical help during measure- 
ments. 

References 



D. C. Johnston, Advances in Physics 59, 803 (2010). 
G. R. Stewart, Rev. Mod. Phys. 83, 1589 (2011). 

M. Daeumling, J. M. Seuntjens, and D. C. Larbalestier, Nature 346, 332 (1990). 
L. Krusin-Elbaum, L. Civale, V. M. Vinokur, and F. Holtzberg, Phys. Rev. Lett. 69, 2280 (1992). 

D. Giller, A. Shaulov, R. Prozorov, Y. Abulafia, Y. Wolfus, L. Burlachkov, Y. Yeshurun, E. Zeldov, V. M. Vinokur, J. L. 
Peng, and R. L. Greene, Phys. Rev. Lett. 79, 2542 (1997). 

D. Giller, A. Shaulov, Y. Yeshurun and J. Giapintzakis, Phys. Rev. B 60, 106 (1999). 

K. Deligiannis, P. A. J. de Groot, M. Oussena, S. Pinfold, R. Langan, R. Gagnon, and L. Taillefer, Phys. Rev. Lett. 79, 
2121 (1997); T. Nishizaki, T. Naito and N. Kobayashi, ibid. 58, 11169 (1998). 

B. Khaykovich, E. Zeldov, D. Majer, T. W. Li, P. H. Kes and M. Konczykowski, Phys. Rev. Lett. 76, 2555 (1996). 
Y. Abulafia, A. Shaulov, Y. Wolfus, R. Prozorov, L. Burlachkov, Y. Yeshurun, D. Majer, E. Zeldov, H. Wiihl, V. B. 
Geshkenbein, V. M. Vinokur, Phys. Rev. Lett. 77, 1596 (1996). 
S. Bhattacharya, M. J. Higgins, Phys. Rev. Lett. 70, 2617 (1993). 

T. Klein, R. Marlaud, C. Marcenat, H. Cercellier, M. Konczykowski, C. J. van der Beek, V. Mosser, H. S. Lee, S. I. Lee, 
Phys. Rev. Lett. 105, 047001 (2010). 

T. Klein, I. Joumard, S. Blanchard, J. Marcus, R. Cubitt, T. Giamarchi, and P. Le Doussal, Nature 413, 404 (2001). 
B. Shen, P. Cheng, Z. Wang, L. Fang, C. Ren, L. Shan, and H. H. Wen, Phys. Rev. B 81, 014503 (2010). 
R. Prozorov, N. Ni, M. A. Tanatar, V. G Kogan, R. T. Gordon, C. Martin, E. C. Blomberg, P. Prommapan, J. Q. Yan, 
S. L. Bud'ko, and P. C. Canfield, Phys. Rev. B 78, 224506 (2008). 

N. Haberkorn, M. Miura, B. Maiorov, G F. Chen, W. Yu, and L. Civale, Phys. Rev. B 84, 094522 (2011). 
S. Salem-Sugui, Jr., L. Ghivelder, A. D. Alvarenga, L. F. Cohen, K. A. Yates, K. Morrison, J. L. Pimentel, Jr., H. Luo, Z. 
Wang, and H. H. Wen, Phys. Rev. B 82, 054513 (2010). 
[17] M. R. Eskildsen, L. Ya Vinnikov, I. S. Veshchunov, T. M. Artemova, T. D. Blasius, J. M. Densmore, C. D. Dewhurst, N. 
Ni, A. Kreyssig, S. L. Bud'ko, P. C. Canfield, A. I. Goldman, Physica C 469, 529 (2009). 

10 



[1] 
[2] 
[3] 
[4] 
[5 

[6] 
[7] 

[8 
[9 

[10 

[11 

[12 
[13 
[14 

[15 
[16 



[18] D. S. Inosov, T. Shapoval, V. Neu, U. Wolff, J. S. White, S. Haindl, J. T. Park, D. L. Sun, C. T. Lin, E. M. Forgan, M. 

S. Viazovska, J. H. Kim, M. Laver, K. Ncnkov, O. Khvostikova, S. Kiihnemann, and V. Hinkov, Phys. Rev.B 81, 014513 

(2010). 
[19] L. Ya Vinnikov, T. M. Artemova, I. S. Veshchunov, N. D. Zhigadlo, J. Karpinski, P. Popovich, D. L. Sun, C. T. Lin, A. 

V. Boris, JETP letters 90, 299 (2009). 
[20] R. Kopeliansky, A. Shaulov, B. Ya Shapiro, Y. Yeshurun, B. Rosenstein, J. J. Tu, L. J. Li, G. H. Cao, and Z. A. Xu, Phys. 

Rev. B 81, 092504 (2010). 
[21] A. K. Pramanik, L. Harnagea, C. Nacke, A. U. B. Wolter, S. Wurmehl, V. Kataev, and B. Biichner, Phys. Rev. B 83, 

094502 (2010). 
[22] S. Salem-Sugui Jr., L. Ghivelder, A. D. Alvarenga, L. F. Cohen, L. Huiqian, and L. Xingye, Phys. Rev. B 84, 052510 

(2011). 
[23] A. K. Pramanik, L. Harnagea, S. Singh, S. Aswartham, G. Bchr, S. Wurmehl, C. Hess, R. Klingeler and B. Biichner, Phys. 

Rev. B 82, 014503 (2011). 
[24] H. Yang, C. Ren, L. Shan, and H. H. Wen, Phys. Rev. B 78, 092504 (2008). 

[25] R. Prozorov, M. E. Tillman, E. D. Mun, and P. C. Canfield, New Journal of Physics 11, 035004 (2009). 
[26] J. D. Moore, L. F. Cohen, Y. Yeshurun, A. D. Caplin, K. Morrison, K. A. Yates, C. M. McGilvery, J. M. Perkins, D. W. 

McComb, C. Trautmann, Z. A. Ren, J. Yang, W. Lu, X. L. Dong, and Z. X. Zhao, Supercond. Sci. Technol. 22, 125023 

(2009). 
[27] C. J. van der Beek, G. Rizza, M. Konczykowski, P. Fertey, I. Monnet, T. Klein, R. Okazaki, M. Ishikado, H. Kito, A. 

Iyo, H. Eisaki, S. Shamoto, M. E. Tillman, S. L. Budko, P. C. Canfield, T. Shibauchi, and Y. Matsuda, Phys. Rev. B 81, 

174517 (2010). 
[28] H. G. Schnack, R. Griessen, J. G. Lensink, and H. H. Wen, Phys. Rev. B 48, 13178 (1993). 
[29] H. H. Wen, H. G. Schnack, R. Griessen, B. Dam, and J. Rector, Physica C 241, 353 (1995). 
[30] D. Bhoi, P. Mandal, and P. Choudhury, Supercond. Sci. Technol. 21 125021 (2008); D. Bhoi, L. S. Sharath Chandra, P. 

Choudhury, V. Ganesan, and P. Mandal, Supercond. Sci. Technol. 22 095015 (2009). 
[31] D. Bhoi, P. Mandal, P. Choudhury, S. Dash, and A. Banerjee, Physica C 471, 258 (2011). 
[32] T. H. Johansen, M. Baziljevich, D. V. Shantsev, P. E. Goa, Y. M. Galperin, W. N. Kang, H. J. Kim, E. M. Choi, M.-S. 

Kim, and S. I. Lee, Europhys. Lett. 59, 599 (2002). 
[33] R. J. Claim, Physica C 50, 153 (1988). 

[34] K. H. Muller, C. Andrikidis, H. K. Liu, and S. X. Dou, Phys. Rev. B 50, 10218 (1994). 
[35] Y. L. Chen, Y. J. Cui, Y. Yang, Y. Zhang, L. Wang, C.H. Cheng, C. Sorrell, and Y. Zhao, Supercond. Sci. Technol. 21, 

115014 (2008). 
[36] A. Yamamoto, A. A. Polyanskii, J. Jiang, F. Kametani, C. Tarantini, F. Hunte, J. Jaroszynski, E. E. Hellstrom, P. J. Lee, 

A. Gurevich, D. C. Larbalestier, Z. A. Ren, J. Yang, X. L. Dong, W. Lu and Z. X. Zhao, Supercond. Sci. Technol. 21, 

095008 (2008) 
[37] Y. Yeshurun, A. P. Malozemoff, and A. Shaulov, Rev. Mod. Phys. 68, 911 (1996). 

[38] H. H. Wen, S. L. Li, Z. W. Zhao, H. Jin, Y. M. Ni, Z. A. Ren, G. C. Che, and Z. X. Zhao, Physica C 363, 170 (2001). 
[39] P. W. Anderson, Phys. Rev. Lett. 9, 309 (1962). 
[40] A. P. Malozemoff, Physica C 185-189, 264 (1991). 

[41] M. V. Feigelman, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Phys. Rev. Lett. 63, 2303 (1989). 
[42] E. Zeldov, N. M. Amer, G. Koren, A. Gupta, R. J. Gambino, and M. W. McElfresh, Phys. Rev. Lett. 62 3093 (1989). 
[43] H. H. Wen, Z. X. Zhao, S. L. Yan, L. Fang, and M. He, Physica C 312, 274 (1999). 

[44] H. Jin, H. H. Wen, H. P. Yang, Z. Y. Liu, Z. A. Ren, G. C. Che, and Z. X. Zhao, Appl. Phys. Lett. 83, 2626 (2003). 
[45] G. Blatter, M. V. Feigelman, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Rev. Mod. Phys. 66, 1125 (1994). 
[46] R. Greissen, H. H. Wen, A. J. J. van Dalen, B. Dam, J. Rector, H. G. Schnack, S. Libbrecht, E. Osquiguil, and Y. 

Bruynseraede, Phys. Rev. Lett. 72, 1910 (1994). 
[47] C. J. van der Beek, M. Konczykowski, S. Kasahara, T. Terashima, R. Okazaki, T. Shibauchi, and Y. Matsuda, Phys. Rev. 

Lett. 105, 267002 (2010). 
[48] G. Blatter, V. B. Geshkenbein, J. A. G. Koopmann, Phys. Rev. Lett. 92 067009 (2004). 
[49] T. Giamarchi, P. Le Doussal, Phys. Rev. B 55, 6577 (1997). 

[50] J. Kierfeld and V. Vinokur, Phys. Rev. B 61, R14928 (2002); J. Kierfeld, V. Vinokur, Phys. Rev. B 69, 024501 (2004). 
[51] G. P. Mikitik, E. H. Brandt, Phys. Rev. B 64 184514 (2001); G. P. Mikitik, E. H. Brandt, Phys. Rev. B 68 054509 (2003); 

G. P. Mikitik, E. H. Brandt, Phys. Rev. B 71 012510 (2005). 
[52] B. Rosenstein, B. Y. Shapiro, I. Shapiro, Y. Bruckental, A. Shaulov, and Y. Yeshurun, Phys. Rev. B 72, 144512 (2005). 
[53] R. Gilardi, J. Mesot, A. Drew, U. Divakar, S. L. Lee, E. M. Forgan, O. Zaharko, K. Conder, V. K. Aswal, C. D. Dewhurst, 

R. Cubitt, N. Momono, and M. Oda, Phys. Rev. Lett. 88, 217003 (2002). 
[54] V. Vinokur, B. Khaykovich, E. Zeldov, M. Konczykowski, R. A. Doyle, and P. Kes, Physica C 295, 209 (1998). 
[55] T. Nishizaki, T. Naito, and N. Kobayashi, Phys. Rev. B 58, 11169 (1998). 

[56] Y. Radzyner, A. Shaulov, Y. Yeshurun, I. Felner, K. Kishio, and J. Shimoyama, Phys. Rev. B 65, 214525 (2002). 
[57] L. Miu, Y. Tanabe, T. Adachi, Y. Koike, D. Miu, G. Jakob, and H. Adrian, Phys. Rev. B 78 024520 (2008); L. Miu, T. 

Adachi, K. Omori, and Y. Koike, D. Miu, Phys. Rev. B 82 064520 (2010). 
[58] D. Miu, T. Noji, T. Adachi, Y. Koike, and L. Miu, Supercond. Sci. Technol. 25 115009 (2012). 
[59] K. Hashimoto, T. Shibauchi, T. Kato, K. Ikada, R. Okazaki, H. Shishido, M. Ishikado, H. Kito, A. Iyo, H. Eisaki, S. 

Shamoto, and Y. Matsuda, Phys. Rev. Lett. 102, 017002 (2009). 

11 



